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ABSTRACT 

This paper addresses the question: how should ./V n-dimensional subspaces of m-dimensional 
Euclidean space be arranged so that they are as far apart as possible? The results of extensive 
computations for modest values of N, n, m are described, as well as a reformulation of the 
problem that was suggested by these computations. The reformulation gives a way to describe 
n-dimensional subspaces of m-space as points on a sphere in dimension (m — l)(m + 2)/2, 
which provides a (usually) lower-dimensional representation than the Pliicker embedding, and 
leads to a proof that many of the new packings are optimal. The results have applications to 
the graphical display of multi-dimensional data via Asimov's "Grand Tour" method. 



1. Introduction 



Although there is a considerable literature dealing with Grassmannian spaces (see for ex- 



ample H], |29|, |45j], |2J, |20|], |18Q, |4J]), the problem of finding the best packings in such 
spaces seems to have received little attention. 

We have made extensive computations on this problem, and have found a number of puta- 
tively optimal packings. These computations have led us to conclude that the best definition 
of distance on Grassmannian space is the "chordal distance" defined in Section 2. 

Sections 3, 4, 5, 6 discuss the problems of packing lines in R 3 , planes in M 4 , n-spaces in 
M m , and lines in R m , respectively. 

The complete results have been placed in a database on NJAS's home page. Our search 
has concentrated on packings of N < 55 subspaces of G(m,n), for m < 16, n < 3. 

2. Grassmannian space 

The Grassmannian space G(m, n) is the set of all n-dimensional subspaces of real Euclidean 
m-dimensional space R m . This is a homogeneous space isomorphic to 0{m) / (O(n) x 0(m—n)), 
and forms a compact Riemannian manifold of dimension n(m — n). 

We first discuss how to define the distance between two n-planes P, Q £ G(m,n). The 
principal angles 0\, . . . ,9 n G [0, n/2] between P and Q are defined by (we follow [p^tl , p. 584) 

cos 9i = max max u ■ v = ui ■ V{ , 

u&P v&Q 

for i = 1, . . . , n, subject to u • u = v ■ v = 1, u ■ Uj = 0, v ■ vj = (1 < j < i — 1). The vectors 
{ui} and {vj} are principal vectors corresponding to the pair P and Q. 

Wong Ej| shows that the geodesic distance on G(m, n) between P and Q is 



d g (P,Q) = ^9{ + -.. + el . (2.1) 

However, this definition has one drawback: it is not everywhere differentiable. Consider the 
case n = 1, for example, and hold one line P fixed while rotating another line Q (both passing 
through the origin). As the angle 4> between P and Q increases from to it, the principal angle 
Q\ increases from to tt/2 and then falls to 0, and is non-differentiable at tt/2 (see Fig. [I]). 

Although one might expect this non-differentiability to be a mere technicality, it does in 
fact cause considerable difficulties for our optimizer, especially in higher dimensions in cases 
when many distances fall in the neighborhood of singular points of d g . 







7T/2 



Figure 1: Principal angle 6\ between two lines as the angle between them increases from to 

7T. 

An alternative measure of distance, which we call the chordal distance, is given by 

d c (P, Q) = ^sin 2 #i + --- + sin 2 6 n . (2.2) 

The reason for the name will be revealed later. This approximates the geodesic distance when 
the planes are close, has the property that its square is differentiable everywhere, and, as we 
shall attempt to demonstrate, has a number of other desirable features. 

A third definition has been used by Asimov || and Golub and Van Loan 
namely 

d m (P,Q) = n = max 6»j . 

i=l,...,n 

This shares the vices of the geodesic distance. 

Of course for n = 1 all three definitions are equivalent, in the sense that they lead to the 
same optimal packings. 

We can now state the packing problem: given N, n, m, find a set of n-planes Pi, ... , P/v £ 
G(m,n) so that mmd(Pi, Pj) is as large as possible, where d is either geodesic or chordal dis- 
tance. Since G{m, n) is compact, the problem is well-defined. Because G{m, n) and G(m, m—n) 
are essentially the same space, we will usually assume that n < m/2. 

We also need some further terminology. A generator matrix for an n-plane P £ G(m, n) is 
an n x m matrix whose rows span P. The orthogonal group 0(m) acts on G(m,n) by right 
multiplication of generator matrices. The automorphism group of a subset {Pi, ■ ■ ■ ,Pn} C 
G(m,n) is the subset of 0(m) which fixes or permutes these planes. 

By applying a suitable element of 0(m) and choosing appropriate basis vectors for the 
planes, any given pair of n-planes P, Q with n < m/2 can be assumed to have generator 



19], p. 584, 
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matrices 
and 

" cos 9\ 


. 

respectively, where 0\, . . . ,0 n are the principal angles between them (|45H, Theorem 2). 
3. Packing lines through the origin in 3-space 

Our initial work on this problem was prompted by a question raised in 1992 by Julian 
Rosenman, an oncologist at the Univ. of North Carolina School of Medicine and Computer 
Science, in connection with the treatment of tumors using high energy laser beams J3f|. He 
asked for the best way to separate N lines through a given point in R 3 , or in other words for 
the best packings in G(3, 1). 

Together with W. D. Smith, we had been carrying out an extensive search for the best 
packings of a given number of points on S 2 , i.e. spherical codes [^f| , fl25| , and we therefore 
modified our programs to search instead for packings of lines. We omit the details of this 
search, since we later recalculated these results using the more general methods described in 
Section 4. 

The results are summarized in Table |l|, which gives, for each value of N in the range 2 to 
55, the minimal angle of the best packing we have found of N lines, and for comparison the 
minimal angle of the best packing known of 2N points on S 2 (taken from We see that 

requiring a packing of 2iV points on S 2 to be antipodal is a definite handicap: only in the cases 
of 6 and 12 points do the antipodal and unrestricted packings coincide. Decimals in the tables 
have been rounded to four decimal places. 

The last two columns of the table specify the largest automorphism groupQ we have found 

of any such best antipodal packing of 2N points. The fourth column gives the order of the 

group and the fifth and sixth columns its name in the orbifold notation (cf. |l0|| ) and as the 

double cover of a rotation group. The symbol ±G indicates that the group consists of the 

matrices ±M for M £ Q, where Q is a cyclic (C), dihedral (2?), tetrahedral (T), octahedral 

(O) or icosahedral (I) group. In each case the subscript gives the order of the rotation group. 
lr That is, the subgroup of 0(3) that fixes or permutes the 2N points. 
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Table 1: Best packings found of N lines through origin in R 3 . (The third column gives the 
best packing known of 2N points on a sphere.) 



AT „ £ 

No. 01 
lines 

AT 

iv 


Mm. 
angle 


Mm. 
angle 
(packing) 


Group 
order 


Group 
name 


Group 
structure 


Notes 


2 


90.0000 


109.4712 


16 


*224 


±T> 8 


square 


3 


90.0000 


90.0000 


48 


*432 


±C 2 4 


octahedron 


4 


70.5288 


74.8585 


48 


*432 


1 /7~\ 

±e> 24 


cube 


5 


63.4349 


66.1468 


20 


2*5 


±£>10 


pentagonal antipnsm (see note) 


6 


63.4349 


63.4349 


120 


*532 


±^60 


icosahedron 


7 


54.7356 


55.6706 


48 


*432 


1 /7~\ 

±C 2 4 


l 1*11 11 

rhombic dodecahedron 


8 


A r\ r» o r\r\ 

49.6399 


f £~\ £~\ A A A 

52.2444 


2 


X 




see note 


9 


47.9821 


49.5567 


6 


3x 




see note 


9 


47.9821 


49.5567 


4 


2* 


1 /I 




10 


46.6746 


47.4310 


24 


*226 


±Pl2 


hexakis bi-antiprism (see note) 


11 


44.4031 


44.7402 


20 


2*5 


±£>10 




12 


41.8820 


43.6908 


48 


*432 


±024 


rhombicuboctahedron 


13 


on a i o. 1 

39.8131 


41.0377 


4 


2* 






14 


38.6824 


39.3551 


2 


X 


±Ci 




15 


38.1349 


38.5971 


20 


2*5 


±£>10 


see note 


16 


37.3774 


37.4752 


120 


*532 


±^60 


j_ 1 ■ 11 1 1 

pentakis dodecahedron 


17 


35.2353 


35.8078 


2 


X 


±Ci 




18 


34.4088 


35.1897 


6 


3x 


±c 3 




19 


33.2115 


34.2507 


2 


X 






20 


32.7071 


33.1584 


8 


*222 


±£>4 




21 


32.2161 


32.5064 


10 


5x 






22 


31.8963 


31.9834 


12 


2*3 


±p 6 




23 


30.5062 


30.9592 


2 


X 


±Ci 




24 


30.1628 


30.7628 


24 


3*2 






25 


29.2486 


29.7530 


6 


3x 






26 


28.7126 


29.1948 


4 


2* 


±c 2 




27 


28.2495 


28.7169 


2 


X 
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Table 1 (cont.): Best packings found of N lines through origin in R 3 . 



No. of 


Min. 


Min. 


Group 


Group 


Groi 


lines 


angle 


angle 


order 


name 


struct 


N 




(packing) 








28 


27.8473 


28.1480 


2 


x 




29 


27.5244 


27 5564 


28 


2*7 




30 


26 9983 


27 1928 


20 


2*5 


- L - i> TU 


31 


26 4987 


26 6840 


10 


5x 


±C^ 


32 


25 9497 


26 2350 


4 


2* 




33 


25 5748 


25.9474 


24 


3*2 


- J -- L LZ 


34 


25 2567 


25 4638 


6 


3x 


±C^ 


35 


24 8702 


25 1709 


2 


x 




36 


24 5758 


24 9265 


6 


3x 


±Ci 


37 


24.2859 


24.4209 


2 


x 




38 


24 0886 


24 1282 

^ J_ > _l_ ^ V... 1 


2 


x 


zbCi 


39 


23 8433 


23 9310 


6 


3x 




40 


23 3293 


23 5531 


2 


x 


=bCi 


41 


22 9915 


23 1946 


2 


x 


=bCi 


42 


22 7075 


23 0517 

a V 1 KJ L 1 


6 


3x 


- LL -'0 


43 


22 5383 

J— • KJ KJ KJ 


22.6744 


12 


2*3 




44 


22.2012 


22.4679 


2 


x 




45 


22.0481 


22.1540 


4 


2* 


±c 2 


46 


21.8426 


22.0276 


2 


X 




47 


21.6609 


21.7221 


2 


X 




48 


21.4663 


21.5206 


6 


3x 


±c 3 


49 


21.1610 


21.3711 


2 


x 


±Ci 


50 


20.8922 


21.0312 


2 


X 


±Ci 


51 


20.6903 


20.8556 


2 


X 


±Ci 


52 


20.4914 


20.6566 


2 


X 


±Ci 


53 


20.2685 


20.4394 


2 


X 


±Ci 


54 


20.1555 


20.3044 


6 


3x 


±c 3 


55 


20.1034 


20.1110 


120 


*532 


±^60 



Notes 



see note 



In some cases the best packings can be obtained by taking the diameters of a known 
polyhedron, and if so this is indicated in the final column of the table. Other entries in this 
column refer to the brief descriptions given below. 

The entries for N < 6 were shown to be optimal by Fejes Toth in 1965 |ll| (see also 
Rosenfeld |p7| ), and the 7- line arrangement will be proved optimal in Section 5. The solutions 
for N > 8 are the best found with over 15000 random starts with our optimizer. There is no 



5 



guarantee that these are optimal, but experience with similar problems suggests that they will 
be hard to beat and in any case will be not far from optimal. 

For N = 1,2,3,6 the solutions are known to be unique, for N = 4 there are precisely 
two solutions (p6|, 36], and for N = 5,7,8 the solutions appear to be unique. For 



larger values of N, however, the solutions are often not unique. For N = 9 lines there are 
two different solutions, and in the range N < 30 the solutions for 10, 22, 25, 27, 29 lines (and 
possibly others) contain lines that "rattle", that is, lines which can be moved freely over a 
small range of angles without affecting the minimal angle. 

Notes on Table 1 

N = 5. Five of the six diameters of a regular icosahedron. 

N = 8. The putatively optimal arrangement forms an unpleasant-looking configuration of 16 
antipodal points on S 2 , with no further symmetry. The convex hull is shown in Fig. 0. 
In contrast, the putatively best packing of 16 points (Fig. 0) has a group of order 16. 



Figure 2: Best antipodal packing found of 16 points. 



N = 9. There appear to be two inequivalent solutions. The nicest has symmetry group 3x, 
of order 12 (or [2 + , 6 + ] in Coxeter's notation [13|) — see Fig. ||. The points lie in six horizontal 
layers of equilateral triangles. The nine points in the Northern hemisphere can be located by 
drawing seven equilateral spherical triangles of edge length 47.98213264° , arranged as in Fig. pi 
The North pole is located at the midpoint of the central triangle. The second solution, shown 
in Fig. |(| has group 2*, of order 4. 
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Figure 3: Best (unrestricted) packing of 16 points known. 




Figure 4: Antipodal packing of 18 points with group of order 6. 

N = 10. This is a "hexakis bi-antiprism" , since it consists of a bi-(hexagonal antiprism) 
together with an axial line that can "rattle" , giving infinitely many solutions. The axis appears 
horizontally in Fig. [?]. 

N = 15. Combinatorially, this is an "axially depleted pentakis dodecahedron", obtained 
from a pentakis dodecahedron (the solution for N = 16) by omitting two opposite vertices (see 
Fig. H). However, the angular separation is slightly greater than for iV = 16. 

N = 55. The 110 antipodal points can be taken to be the union of the vertex sets of a 
dodecahedron (20), an icosidodecahedron (30) and a truncated icosahedron (60) — see Fig. ^. 
The 15 lines through the icosidodecahedral points "rattle". This arrangement of points can be 
obtained from a "geodesic dome" or 12-fold reticulated icosahedron |l^| by omitting its twelve 
pentagonal vertices, and so may be described as consisting of the vertices of a "depleted 12-fold 
reticulated icosahedron". Incidentally, our best solution for 61 lines has a symmetry group of 
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Figure 5: 




Figure 6: Antipodal packing of 18 points with group of order 4. 

order 6, and is not found by placing points at the centers of the pentagonal faces of Fig. ^. 
Packings of lines in higher dimensional spaces will be discussed in Section 6. 

4. Packing two-dimensional planes in R 4 

In H Asimov proposed a technique called the "Grand Tour" for displaying multi-dimensional 
data on a two-dimensional computer screen. His idea was to choose a finite sequence or "tour" 
of two-dimensional planes that are in some sense "representative" of all planes, and to project 
the data onto each plane in turn, in the hope that the viewer will be able to notice any pattern 
or structure that is present. This technique has been implemented in the XGobi program ||, 

@- 

In 1993 Dianne Cook asked us if we could modify our algorithm for finding spherical codes, 
in order to search for packings in G(m, 2) for m > 4. Furthermore, for the Grand Tour 
application, there should be a Hamiltonian circuit through the planes, in which edges indicate 
neighboring planes. We accomplished this in the following way. 
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Figure 7: Best antipodal packing found of 20 points. 



Figure 8: Best antipodal packing found of 30 points. 

Following the methods we had used to find packings in S d and experimental designs in var- 
ious spaces (see p2f , Sec. 4.11; p3~| , pip , corresponding to a set of planes S = {Pi, . . . , P/v} C 
G(m, n) we define a potential 

*° (a) -§ *w,q)-A - 

where ^4 is a suitably chosen constant. There is a similar definition for § g involving d g . Initially 
A is set to and S is a randomly chosen set of planes. We invoke the Hooke- Jeeves "pattern 
search" optimization algorithm to modify S, attempting to minimize $> C (S). After a fixed time 
(we used 100 steps of the optimizer), A is advanced halfway to the current minimal distance 
between the Pi, and the process is repeated, terminating when no further improvement is 
obtained to the accuracy of the machine. The whole process is repeated with several thousand 
random starts — and also with initial configurations taken from other sources, such as subsets 
or supersets of other arrangements — and the best final configuration is recorded. 

The partial derivatives of & C (S) were found analytically, but for $> g (S) we calculated them 
by numerical differentiation. 
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Figure 9: Best antipodal packing found of 110 points. 



We began by computing a table of packings of N planes in G(4, 2), for N < 50, with 
respect to both d c and d g . The results for packings of N < 24 planes are summarized in 
Table ||. Columns 2 and 3 refer to the best packings found for the chordal distance, and give 
both 

d\ = min d\{P h Pj ) and dl = min dl {P { , Pj ) 

for such packings. Columns 4 and 5 refer similarly to the best packings found with respect to 
geodesic distance. 

Understanding these results was hindered by the fact that planes are, as the name suggests, 
plain, with no distinguishing features, and when produced by the computer appear as random 
generator matrices referred to a random coordinate frame. However, we found that the set 
of principal vectors in the planes could often be used to find a coordinate system that would 
reveal some of the structure of the planes. 

For two planes the best packing is the same for both metrics, and consists of two mutually 
orthogonal planes with principal angles vr/2, 7r/2, so S c = 2, dg = 7r 2 /2. 

For three planes the two answers are different. The best packing for chordal distance 
(Fig. |i~C|a) has principal angles 7r/3, vr/3 between each pair of planes, so d\ = 3/2, d 2 g = 2vr 2 /9, 
whereas the best packing for geodesic distance (Fig. [Tofb) has angles 7r/6, 7r/2, so d 2 c = 5/4, 
d 2 g = 5vr 2 /18. 

Postponing discussion ol N = A and 5 for the moment, let us consider the case of six 
planes, where we discovered that the answer for both metrics formed a regular simplex^: the 

principal angles between every pair of planes were arcsin l/\/5 and ir/2. This set of planes is 

2 Although not an isoclinic configuration, cf. Wong Q|. There seems to be only a slight overlap connection 
between that problem and ours. 
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Table 2: Comparison of putatively best packings of N planes in G(4, 2) with respect to d c and 

dg. 



Best w.r.t. d c Best w.r.t. d g 



N 




dl 


4 
y 




dl 




4 

y 


2 


2 


.0000 


4.9348 


2. 


,0000 


4. 


,9348 


3 


1 


.5000 


2.1932 


1. 


,2500 


2. 


.7416 


4 


1 


,3333 


1.8253 


1 


,2000 


2. 


.6824 


5 


1 


.2500 


1.9739 


1 


,2000 


2 


.6824 


6 


1 


,2000 


2.6824 


1 


,2000 


2. 


.6824 


7 


1 


,1667 


1.6440 


0. 


.9875 


2 


,1281 


8 


1 


.1429 


1.5818 


0. 


.9700 


1 


.9235 


9 


1 


.1231 


1.5175 


o. 


.9501 


1 


.8087 


10 


1 


.1111 


1.5725 


0. 


.9764 


1 


.7886 


11 


1 


.0000 


1.2715 


0. 


,9247 


1 


,6711 


12 


1 


.0000 


1.3413 


0. 


.9204 


1 


.6416 


13 


1 


,0000 


1.2348 


0. 


.9133 


1 


.5600 


14 


1 


.0000 


1.2337 


o. 


.8933 


1 


.5327 


15 


1 


,0000 


1.2337 


0. 


.8923 


1 


,5284 


16 


1 


.0000 


1.2337 


0. 


.8904 


1 


,5210 


17 


1 


,0000 


1.2337 


0. 


.8549 


1 


,3925 


18 


1 


.0000 


1.2337 


0. 


.8504 


1 


.3768 


19 


0. 


.9091 


1.1666 


0. 


.8412 


1 


.3477 


20 





.9091 


1.1666 


0. 


.8351 


1 


,3284 


21 





.8684 


1.0352 


0. 


.8225 


1 


.2834 


22 


0. 


.8629 


1.0592 


0. 


.8046 


1 


,2385 


23 


0. 


.8451 


1.0081 


0. 


.7910 


1 


,2012 


24 


0. 


.8372 


0.9901 


0. 


.7812 


1 


.1707 



conveniently described using simplicial coordinates. Let A, B, C, D, E be the vectors from 
the center of a regular simplex in R 4 to its vertices, with A + B + C + D + E = 0, and write 
[a, b, c, d, e] for a A + bB + cC + dD + eE. Then one of the six planes is spanned by [1, r, 1, 0, 0] 
and its cyclic shifts, where r = (1 + \/5)/2, and the other five planes are obtained from it by 
taking even permutations of these coordinates. (We later found an equally good packing with 
respect to chordal distance that was not a simplex with respect to geodesic distance. This is 
described below.) 

The six planes intersect the surface of the unit ball of R 4 in a remarkable link, the 
"Hexalink" . It can be shown that apart from the n-component unlinks and Hopf links, which 
exist for all n, the Hexalink is the only link of circular rings in which there are orientation- 
preserving symmetries taking any two links to any other two. (Details will be given elsewhere.) 
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1 





1 








1 





1 


1 





1 

2 


73 
2 





1 


V3 
2 


1 
2 


1 





1 

2 


73 
2 





1 


73 
2 


1 
2 


(a) 



1 


o 


o 


o 





1 

72 





1 

72 





1 














1 

V2 


1 








1 





1 

75 








1 

72 



(b) 

Figure 10: Generator matrices for best packings of three planes in G(4, 2) with respect to (a) 
chordal and (b) geodesic distance. 

The discovery of this 6-vertex simplex was initially somewhat of a surprise, since G(4, 2) 
is only a 4-dimensional manifold. It did suggest that the answer should somehow be related 
to the six equi-angular diameters of the icosahedron, and led to the following reformulation of 
the problem. 

We remind the reader that any element a of SO (4) may be represented as 

a : x h 



where x = XQ + x\i + X2j + x$k represents a point on S s and £, r are unit quaternions [15]. The 
pair —£, —r represent the same a. The correspondence between a and ±(£,r) is one-to-one. 

Given a plane P E £7(4, 2), let a be the element of SO (4) that fixes P and negates the points 
of the orthogonal plane P ± . Then a 2 = 1, and for this a, it is easy to see that £ = £ii+£2j+£sk 
and r = r\i + r<i3 + r^k are purely imaginary unit quaternions. This establishes the following 
result (which can be found for example in Leichtweiss |29f] ). 

Theorem 1. A plane P G G(4, 2) is represented by a pair (£,r) G S 2 x S 2 , with (—£,—r) 
representing the same P. 

There are simple formulae relating P and (£, r), pointed out to us by Simon Kochen. Given 
(£, r), if £ 7^ — r then P is spanned by the vectors corresponding to the quaternions u = 1 — £r 
and v = £ + r. The special case when £ = r is even simpler: take u = 1, v = £. If £ = —r, take 
u and i> to be purely imaginary unit quaternions such that £, u, v correspond to a coordinate 
frame. Conversely, if P is spanned by two orthogonal unit vectors represented by quaternions 
u, v, then ±(£, r) = ±(uv — vu, vu — uv). 

Given two planes P,Q £ G(4, 2), represented by ±(£, r), ±(£' , r'), respectively, the principal 
angles #i, 62 between them may be found as follows. Let (ft (resp. tp) be the angle between £ 

12 



and £' (resp. r and r'), with < eft, if) < it. If <p + "0 > ft, replace (f> by n — (f> and if) by 7r — if), 
so that < <j) + ?/> < 7T, with (j> < if) (say). Then 



(4.1) 



d 2 9 (P,Q) 
d 2 c (P,Q) 



^ 2 + 2 



1 — COS if) COS ( 



(4.2) 
(4.3) 



We omit the elementary proof. 

A set S = {Pi, . . . ,Pn} Q G(4,2) is thus represented by a "binocular code" consisting of 
a set of pairs r^) S S 2 x S 2 . We call the list of 2N points ±£i (they need not be distinct) 
the "left code" corresponding to S, and the points ±rj the "right code". Conversely, given 
two multisets L C S 2 , R C S 2 , each of size 2N and closed under negation, and a bijection or 
"matching" / between them that satisfies f(—£) = —f(£), £ E L, we obtain a set of N planes 
in G(4,2). 

The binocular codes for the d c -optimal packings of N = 2, ... ,6 planes are shown in 



Figs. |11|, |12|. Except for N = 3, the left and right codes are identical. Matching points from 
the left and right codes are labeled with the same symbol. For N < 5 the points lie in the 
equatorial plane, and for N < 4 there are repeated points. The points lie on regular figures, 
except for N = 4 where the points are ±(1, 0, 0), ± 



'0,1,2 




N=2 




N=3 




Figure 11: Binocular codes describing best packings of N = 2,..., 5 planes in G(4,2) for 
chordal distance. 
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For N = 6 the left and right codes consist of the 12 vertices of an icosahedron (Fig. |T^). 
Let these be the points 

A(0,±l,±r), A(±r,0,±l), A(±l,±r,0) , (4.4) 



where A = 1/y/r + 2. The matching is obtained by mapping each point to its algebraic 
conjugate (i.e. replacing \/5 by —y/E), and rescaling so the points again lie on a unit sphere. 
As already mentioned, the principal angles between each pair of these planes are arcsin l/\/5 
and tt/2, so d 2 c = 6/5, dl = 2.6824. 





Figure 12: Best packing of 6 planes in G(4, 2) with respect to both metrics. Left and right codes 
comprise vertices of icosahedron. Adjacent vertices in one code are matched with non-adjacent 
vertices in the other code. 

There is a second set of six planes with d 2 = 6/5, but with d 2 only equal to 2.0030. 
Here the left and right codes form the vertices of what we shall call the "anti-icosahedron" , 
consisting of the points A(0,±l,±r), A(±r, 0, ±1), A(±r, ±1,0) (compare fl4.4| )), and shown 
in Fig. [ll|. Topologically this is a "parallel bi-slit cuboctahedron" , obtained by dividing two 
opposite square faces of a cuboctahedron into two triangles by parallel lines. The matching 
sends A(0, ±1, r) to A(0, ±r, 1), A(r, 0, ±1) to A(l, 0, ±r) and A(±r, 1, 0) to A(±l, r, 0). 

Four larger putatively d c -optimal packings that the computer found are also worth men- 
tioning, those of 10, 18, 48 and 50 planes. 

For the 10-plane packing the left and right codes consist of the vertices of a decagonal prism 
with coordinates 

± C ° S r0, Shl r0, VT) ' T = °' ' ' ' ' 9 ' 
where 9 = 7r/5. A typical point (>/2/3 cos r9, y/2/3 sin r9, y/1/3) in the left code is matched 
with (-l) r (y/2/3 cos 3rd, ^/2/3sm3r6, y/ljS) in the right code. Then d\ = 10/9, d 2 = 
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Figure 13: Top and side views of anti-icosahedron. 



1.5725. The most interesting property of this configuration is that, although three different 
sets of canonical angles occur, the chordal distance between every pair of planes is the same 
— this is a regular simplex! 

We find it surprising that this is superior to the packing of ten planes obtained by matching 



the vertices of a dodecahedron to their algebraic conjugates, as shown in Figure |14. 





Figure 14: Matching of vertices of dodecahedron. 



18 planes. Let O denote the set of vertices of a regular octahedron. The binocular code for 
the 18-plane packing is O x O := {(£, r) : £, r 6 O}. Projectively, the 18 two-spaces are the 

edges of a desmic triple^] of tetrahedra (4^] , |[7j] , [|T| . We were pleased to see this configuration 

3 Two tetrahedra are called desmic if they are in perspective from four distinct points — these form a third 
tetrahedron desmic with either. The 18 edges of such a desmic triple are also the edges of a unique other desmic 
triple. 
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appear, since it was already familiar in the context of quantum logic Q. This set of planes can 
be constructed directly by taking the planes spanned by all pairs of mutually perpendicular 
minimal vectors of the D4 lattice (cf. ||]). In this form the planes have generator matrices 
such as 



+ + 00 
+ - 



(12 planes) 



+ + 00 
+ -00 



(6 planes) . 



Three different sets of principal angles occur, 0, n/2; vr/4, 7r/4; and tt/2, tt/2; so that d 2 = 1, 
dg = 1.2337. The automorphism group of this packing has structure [3,4, 3]. 2 and order 2304. 
These 18 planes can also be described as the set of planes that meet the 24-cell in one of its 
equatorial squares. 

48 and 50 planes. Let C denote the set of vertices of a cube. Then the binocular codes for 
the 48- and 50-plane packings are respectively O x C U C x O and OxOuCxC. For both 
arrangements we have d 2 = 2/3, d 2 = 0.7576. Both are believed to be optimal. (We mention 
the 48-plane packing because of its greater symmetry.) 



The matching problem. In order to maximize the minimal chordal distance between the 
planes, the matching between the left and right codes should (from fl4.3| )) be chosen so as to 
minimize the maximal value of cos ip cos 4>. Stated informally, the matching should be such 
that if two points are close together on one of the spheres then the points to which they are 
matched should be far apart on the other sphere. In the case of the icosahedron, for example 
(see Fig. 12), the matching sends adjacent vertices to non-adjacent vertices. There is a unique 
way to do this. 

At this point we were tempted to see if any new record packings in G(4, 2) could be obtained 
by taking the 2N antipodal points corresponding to a good iV-line packing in G(3, 1), and 
matching them with themselves in an optimal way. However, it was not easy to see how to 
solve the matching problem. Fortunately David Applegate (personal communication) found 
that it could be reformulated as an integer programming problem, as follows. 

Given an antipodal set S = {Pi, ■ ■ ■ , P2N} f= S 2 , we wish to find a permutation / of 
S with the property that /(— Pj) = /(Pi), for all i, and such that the minimal value of 
1 — Pi - Pjf(Pi) ■ f(Pj) (i 7^ j) is maximized. If we do the maximization by binary search, we 
may define 7r(Pi,Pj) = 1 if f(Pi) = Pj, or otherwise, with the constraints 

iriPuPj) = TT(-P,-Pj) , 
16 



2N 

J2*(Pi,Pj) = 1, l<i<2N, 

3=1 
2N 

J2ir(Pi,Pj) = 1, l<i<2iV, 
i=i 

and 

ir(Pi, Pj) + ir(P k , Pi) < 1 for all 1 < k,£<2N such that 1 - P { ■ P k Pj ■ P e < M . 

There is a feasible solution if and only if d% > M. 

Applegate kindly implemented this procedure, and used it to solve the matching problem 
for our best packings in G(3, 1) and for various polyhedra. Unfortunately no new records have 
yet been obtained by this method. 

Best packings for the geodesic metric. So far we have mostly discussed packings that 
attempt to maximize the minimal chordal distance between planes. We also computed packings 
for the geodesic distance, and we shall now describe some of them. In general, however, these 
are much less symmetric than the chordal-distance packings, especially for more than 16 planes. 

For N = 4,5 the solutions are subsets of the 6-plane packing. For N = 7 the left code is a 
heptagonal anti-prism and the right code is an equatorial 14-gon. 

For N = 12 the binocular code is the set {±(£,r)}, where i (resp. r) runs through the 
vertices of a regular tetrahedron (resp. equatorial equilateral triangle). 

For N = 16 the left and right codes are the union of two similarly oriented square prisms. 
For larger N the geodesic packings do not seem so interesting. For N = 18, for example, the 
best packing has no non-trivial symmetries. 

Hamiltonian paths. One further question remains to be discussed, that of arranging the 
planes in a circuit in such a way that adjacent planes are close together, for the "Grand 
Tour" application. This turned out to be a much easier problem than finding the packings. 
We handled it in two different ways. For configurations such as the 48-plane arrangement, 
where there was an obvious notion of adjacency — in this case, defining two planes to be 
adjacent if the principal angles are 7r/4, 7r/4 — we represent the packing by a graph with 
nodes representing the planes, and look for a Hamiltonian cycle. In less regular cases, we 
convert the packing into a traveling salesman problem using chordal distance to define the 
distance between nodes, and look for a minimal length circuit. In both cases we were able 
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to make use of the travelling salesman programs of Applegate et al. ||, which can handle 
100-node graphs without difficulty. Some of the packings in the database mentioned in Sect. 1 
(those with suffix .ham) have been arranged in cycles in this way. 



The above reformulation in terms of binocular codes applies only to G(4, 2) (we realized 
from the beginning that this case would be special, since G(4, 2) is the only Grassmannian 
space where the Riemannian metric is not unique |^9|). In the next section we describe a 
second reformulation that applies to the general case. 

5. Packing n-planes in E m 

Three observations contributed to the second reformulation. 

(i) We noticed (see Table 2) that for several examples of N-plane packings in G(4, 2) the 
largest value of d 2 c that we could attain was N/(N — 1), and that in every case this was an 
upper bound. Further experimentation with other packings in G(m, n) for m < 8 led us to 
guess an upper bound of 

, 9 n(m — n) N . 

which again we could achieve for some small values of N. The form of (^^) was suggestive of 
the Rankin bound for spherical codes || or the Plotkin bound for binary codes p3]. 

(ii) Investigation of the 18-plane chordal-distance packing revealed that, with respect to 
chordal distance, this has the structure of a regular orthoplex (or generalized octahedron) with 
18 vertices. Combining this with the fact that the 10-plane packing formed a regular simplex, 
we had strong evidence that (3(4, 2) should have an isometric embedding (with respect to d c ) 
into M 9 . 

(hi) A computer program was therefore written to determine the lowest dimensions into 
which our library of packings in G(m, n) could be isometrically embedded. More precisely, for 
a given set of iV planes in G(m,n), we searched for the smallest dimension D such that there 
are N points in R D whose Euclidean distances coincide with the chordal distances between our 
planes. 

The results were a surprise: it appeared that G(m, n) with chordal distance could be 
isometrically embedded into MP , for D = C™^ 1 ) — 1, independent of n. Furthermore the points 
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representing elements of G(m, n) were observed to lie on a sphere of radius \/n(m — n)/2m in 

Aided by discussions with Colin Mallows, we soon found an explanation: just associate to 
each P G G(m, n) the orthogonal projection map from R m to P. If A is a generator matrix for 
P whose rows are orthogonal unit vectors, then the projection is represented by the matrix 

V = A tr A . (5.2) 

V is an m x m symmetric idempotent matrix, which is independent of the particular orthonor- 
mal generator matrix used to define it. Changing to a different coordinate frame in R m has 
the effect of conjugating V by an element of 0(m). With the help of ( |2.3| ), we see that 

trace V = n . (5-3) 

Thus V lies in a space of dimension C 11 ^" 1 ) — 1. 

Let || || denote the Z/2-norm of a matrix: if M = (M{j), 1 < i,j < m, 



\M\ 



\ 



m m 

Mij = Vtrace M tr M . 

i=l j=i 



For P,Q £ G(m, n), with orthonormal generator matrices A, B, and principal angles 9i, . . . , 6 n , 
an elementary calculation using (|2.3|) , (|2.4| ) shows that 



d 2 c (P, Q) = n - (cos 2 0i -\ h cos 2 



= n - trace A tr AB tr B 

= l\W-Q\\ 2 , (5.4) 

where V, Q are the corresponding projection matrices. 

Note that if we define the "de-traced" matrix V = V — ^I m , then trace V = 0, and 
||"P|| 2 = n ~ ■ We have thus established the following theorem. 

Theorem 2. The representation of n-planes P £ G(m,n) by their projection matrices V gives 
an isometric embedding ofG(m,n) into a sphere of radius y/n(m — n)/n inR D , D = ( m j" ) — 1, 
withd c (P,Q) = -^\\p-Q\\. 

Thus chordal distance between planes is 1/V2 times the straight-line distance between the 
projection matrices (which explains our name for this metric). The geodesic distance between 
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the planes is l/v2 times the geodesic distance between the projection matrices, measured 
along the sphere in M. D . 

Figure |B| attempts to display the embeddings of G(m, 0), G(m, 1), . . . , G(m, m) in R D+1 . 
Since \\V — ^I m \\ 2 = m/4, all the planes lie on the large sphere, centered at \l m , of radius 
y/m/2. The members of G(m, n) lie on the intersection of the large sphere with the plane 
trace('P) = n, which intersection is itself a sphere in H D of radius y/n(m — n)/m centered at 
^I m - A plane P and its orthogonal complement P 1 - are represented by antipodal points on 
the large sphere. 




Figure 15: Embedding of G(m, 0), . . . , G(m, m) into large sphere of radius \/m/2 in Euclidean 
space of dimension m(m + l)/2. G(m,n) lies on sphere of radius \/n(m — n)/m in MP , D = 
(m- l)(m + 2)/2. 

In contrast to this result, we briefly remark, without giving details, that there is no way 
to embed G(m, n) into Euclidean space of any dimension so that the geodesic distance d g on 
G(m, n) is represented by Euclidean distance in that space. Of course the Pliicker embedding, 
in which members of G(m, n) are represented by points in projective space of dimension (™) — 1, 
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also does not give a way to realize either d c or d g as Euclidean distance. (Nor does the Nash 
embedding theorem J3|J|.) Note also that the dimension of the P Kicker embedding is in general 
much larger than the dimension of our embedding. 

Since we have embedded G(m,n) into a sphere of radius y/n(m — n) /m in M. D , we can 
apply the Rankin bounds for spherical codes j35|, and deduce: 



Corollary, (i) The simplex bound: for a packing of N planes in G(m,n), 

l9 n(m — n) N , 

dl < ^ 1 ■ W ~r ■ 5.5 

m N — 1 

Equality requires N < D+l = ( m + 1 ), and occurs if and only if the N points in M> D corresponding 
to the planes form a regular 'equatorial' simplex, 
(ii) The orthoplex bound: for N > ( m ^ _1 ), 

r9 n(m — n) , 

d 2 c < -i L . 5.6) 

m 

Equality requires N < 2D = {m — l)(m + 2), and occurs if the N points form a subset of the 
2D vertices of a regular orthoplex. If N = 2D this condition is also necessary. 



Lemmons and Seidel ( f3l|| , Theorem 3.6) give a bound for equi-isoclinic packings in G(m, n) 
which agrees with (|5.5| ); of course our bound is more general. The case n = 1 of (|5.5|) is given 
in Theorem 3 of |pq ]. 

We were happy to obtain confirmation of (^^). Since d\ can never exceed n, we can also 
write 

o f n(m — n) N 1 , 
^<min^n,— '-— — -\. (5.7) 



m N - 1. 

The Corollary allows us to establish the optimality of many of our packings. In the range 
m < 16, n < 3, iV < 55, there are over 750 cases which appear to meet (|5.7|) or (|5.6| ), in the 
sense that the ratio of d\ to the bound is greater than .9999999 . For n = 2, for instance, the 
following are the packings of N planes in G(m, 2) for m < 10 that appear to achieve either 



fl5.7| ) (listed before the semicolon) or (5^) (after): 



m 


= 4 


2 - 


8, 10; 11 - 


18 


m 


= 5 


4- 


11; 16, 17 




m 


= 6 


3 - 


14; 22, 23 




m 


= 7 


6 - 


18; 29 




m 


= 8 


4- 


21, 28; 37 


-44 


m 


= 9 


6 - 


25; 




m = 


= 10 


5 - 


28; 





21 



However, this is not as meaningful as it at first seems. Consider the packings of N planes 
in G(10, 2). The ratios of d 2 c to the simplex bound for N = 21, ... , 33, rounded to 14 decimal 
places, are as follows: 

21 0.99999999999999 

22 1.00000000000000 

23 0.99999999999999 

24 0.99999999999999 

25 1.00000000000000 

26 0.99999999999999 

27 1.00000000000000 

28 0.99999999921354 

29 0.99999470290071 

30 0.99998661084736 

31 0.99961159256647 

32 0.99909699728979 

33 0.99854979246655 

From this it seems very likely that there exists a 27-plane packing meeting the simplex bound, 
but certainly further investigation is needed to determine if there is a 28-plane simplex packing 
in the neighborhood of the computer's approximate solution. We have carried out such analyses 
in all the putative cases of equality for packing lines (see Sect. 6), for planes in G(4, 2) (as 
already discussed) and for 70 planes in G(8, 4) (see below). The other cases remain to be 
investigated. 

In principle there is no difficulty in settling these questions. All the chordal distances 
between the planes are specified. So we could simply set up a (large) system of quartic equations 
with integer coefficients for the entries in the generator matrices, and ask if there is a real 
solution. 

The chief difficulty in attempting to understand the computer-generated packings is that 
there are usually very large numbers of equally good solutions. The case of 18 planes in G(4, 2) 
is typical in this regard. There is one exceptionally pleasing solution (described in Section 4), 
in which only three different pairs of principal angles occur. However, it seems that there is a 
roughly 29-dimensional manifold of solutions, in which a typical solution, although still having 
d 2 . = 1, has an apparently random set of principal angles. (Our investigation of this question 
is not yet completed.) 

The case of subspaces of dimension n in R 2n or R 2n+1 is especially interesting. The largest 
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possible arrangements of planes that could achieve the two bounds are: 



m n N (simplex) iV(orthoplex) 



2 
3 
4 
5 
6 
7 



3V 

10 V 

15 

21 

28 V 
36 



4V 
10 

18 y/ 
28 
40 
54 

70 y/ 



(5.8) 



Checks indicate that such a packing exists. It is known that the orthoplex bound cannot 
be achieved by 10 planes in G(3, 2), while the other cases are undecided. Our computer 
experiments strongly suggest that no set of 15 planes meets the simplex bound in G(5, 2). 

On the other hand it is possible to find packings of 70 planes in G(8, 4) meeting the bound 
([5.6|). With a considerable amount of effort (a long story, not told here), we determined 
several examples, of which the following is the most symmetrical. Let the coordinates be 
labeled oo, 0, 1, . . . , 6, and take two 4-planes generated by the vectors 



{10000000, 01000000, 00100000, 00001000} 
{11000000, 00101000, 00010001, 00000110} 



(5.9) 



respectively. We obtain 70 planes from these by negating any even number of coordinates, 
and/or applying the permutations (0123456), (oo0)(16)(23)(45) and (124)(365). The principal 
angles are 0,0, f , |; f , f , f , f ; or |, |, |, |, so d? c = 2, d 2 g = 7r 2 /4 (this is not even a local 
optimum with respect to geodesic distance). The full group, which is transitive on the planes, 
is generated by the above operations and by the 8x8 Hadamard matrix 



1 

71 



+--+-++ 
--+-+++ 
-+-+++- 
+-+++ — 
-+++--+ 
+++--+- 
+ + -- + -+. 



(— means —1 
+ means +1) 



It has shape 2 8 _4g and order 5160960 (A n denotes an alternating group of degree n). 

A set of 28 planes in G(7, 3) meeting the simplex bound can be obtained as follows. Label 
the coordinates 0,... ,6, and let v^r have components 1 at 2 r (mod 7), ±\/2 at 3.2 r (mod 7), 
and 0's elsewhere, for r = 0, 1, 2. Then the vectors v^, vf, (product of signs is even) span 
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four planes, and the full set of 28 is found by cycling the seven coordinates^ 

On the other hand, in spite of much effort, we have not been able to find a set of 40 planes 
in G(6, 3) that meet the bound. If this were possible, we would obtain d 2 c = 1.5, whereas 
the best we have been able to achieve is 1.49977. We can get 1.5 with N = 34 planes, and 
1.49977 with 35 < N < 40, but only 1.4297 with 41 planes, which suggests that 1.5 might be 
possible with 40 planes. The discussion given earlier shows that this question could in principle 
be settled by seeing if a certain set of quadratic equations with integer coefficients has a real 
solution. One formulation leads to 500 quartic equations in 360 unknowns, so this approach 
seems hopeless at the present time. 

It is not difficult to invent ways to construct sets of subspaces. Several examples are given 
below. However, many promising ideas have proved useless when confronted with the results 
found by our program. For this reason we have included two extensive tables of the best 
chordal distances that we have found. Table ^ gives d\ for packings in G(m, 2), with m < 10. 
This will serve as a standard against which readers can test their own constructions. 

Table 3 about here 

General constructions. The following are some promising general constructions. 

(i) A skew-symmetric conference matrix^] of order 4a yields a set of 4a unit vectors in C 2a 
with Hermitean inner products ±i/\/4a — 1 ( |14[| , Example 5.8), hence a set of 4a planes in 
G(4a, 2) with d 2 c = 4(2a — l)/(4a — 1) that meet the simplex bound. 

(ii) Use the planes defined by the n-faces of a regular m-dimensional polytope, or the 
Voronoi or Delaunay cells of a lattice, etc. The initial results from this idea have been dis- 
appointing. For example, the 96 two-dimensional faces of the 24-cell in R 4 define 16 different 
planes, forming a packing with d\ = 8/9, inferior to the best such packing (cf. Table |^). 

(iii) Use the minimal vectors in a complex or quaternionic lattice to obtain packings in 
G(2a, 2) or G(4a,4). For example, the 54 minimal vectors of the lattice regarded as a 
three-dimensional lattice over the Eisenstein integers (||, p. 127), produce a packing of nine 

planes in G(6, 2) that meets the simplex bound. 

4 Note added April 1996. P. W. Shor and the third author have recently discovered that the packing of 70 
planes in G(8, 4) can be generalized to a packing of m 2 + m — 2 planes in G(m,m/2) meeting the orthoplex 
bound, whenever m > 2 is a power of 2; and that the packing of 28 planes in G(7, 3) can be generalized to a 
packing of p(p + l)/2 planes in G(p, (p — l)/2) meeting the simplex bound, whenever p is a prime which is either 
3 or congruent to —1 modulo 8. These constructions will be described elsewhere. 

5 For a list of the known orders of such matrices, see Table 7.1 of |40j]. 
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(iv) Restrict the search to generator matrices of O's and l's (as in (|5,9|)) or +1' and — l's, or 
even to rows that are blocks in some combinatorial design, or vectors in some error-correcting 
code (see the example in Section 6). 

(v) Let C be an error-correcting code of length m over GF(2 n ), for example a Reed- 
Solomon code (cf. [^]). If we regard GF(2 n ) as a vector space of dimension n over GF(2), 
each codeword yields annXm matrix whose elements we may take to be +l's and —l's. 
After discarding those of rank less than n, and weeding out duplicates, we obtain a packing in 
G(m,n). The hexacode (||, p. 82), for example, a code of length 6 over GF{4) containing 64 
codewords, produces 28 distinct planes in G(6, 2) with d 2 = 3/4. Unfortunately Table ^| shows 
that the record is 1.2973. 

(vi) Choose a group G with an m-dimensional representation, and a subgroup H of index 
N with an n-dimensional irreducible representation. Find an n-dimensional subspace V C R m 
invariant under H, and take its orbit under G. Many of the conference matrix constructions 
of line-packings described in Section 6 are of this type (using G = L2(q)). 

In some cases we have also used the optimizer to search for packings with a specified 
group. The following is a packing of 28 planes in G(8, 2) meeting the simplex bound that is 
an abstraction of a configuration found by the computer when searching for packings in R 8 
invariant under the permutation (0)(1, 2, 3, 4, 5, 6, 7). Let R : C 4 — > M 8 map (^1,^2,^3,^4) to 
(Reui,Im7;i,Ref2, • • • ,1m. v^), and let a = e 2m / 7 . The 28 planes are spanned by the following 
pairs of vectors: 

R(0,a k ,a 2k ,a Ak ), R{0,ia k ,ia 2k ,ia 4k ), < k < 6 , 
R(l,0,a 2k ,a 4k ), R{i,0,ia 2k ,-ia 4k ), < k < 6 , 
#(1,0^,0,0^), R(i,-ia k ,0,ia 4k ), < k < 6 , 
R(l,a k ,a 2k ,0), R(i,ia k ,-ia 2k ,0), 0<k<6. 

(The pattern of zeros and signs here suggests the Tetracode |J, p. 81.) 

An application: apportioning randomness. The following is a potential application 
of packings in G(m, n) for use in apportioning randomness, in the sense of producing large 
numbers of approximately random points from a few genuinely random numbers. We illustrate 
using the example of 70 planes in G(8,4) described above. Let Ai,...,A?o be generator 
matrices for them. Suppose an exclusive resort wishes to distribute random points in IR 
to its 70 guests, perhaps for use as garage door keys. Let x = (xi, . . . ,xg) be a vector of 
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eight independent Gaussian random variates. Then the hotel would assign y, = AiX tr to its 
i-th guest. By maximizing the chordal distance between the planes we have minimized the 
correlation between the m. 

6. Packing lines in higher dimensions 

Table ^ shows the maximal angular separation found for packings of iV < 50 lines in G(m, 1) 
for m < 9. 

The following packings of N lines in G(m, 1) for m < 10 achieve either the simplex or 
orthoplex bounds: m = 3 : N = 3, 4, 6, 7; m = 4 : N = 4, 5, 11, 12; m = 5 : iV = 5, 6, 10, 16; 
m = 6 : iV = 6,7, 16,22; m = 7 : AT = 7, 8, 14, 28; m = 8 : N = 8, 9; m = 9 : 2V = 
9, 10, 18, 46 - 48, m = 10 : iV = 10, 11, 16. Most of these are described below. 

As mentioned in Section 3, the best packing of 5 lines in G(3, 1) is a subset of the best 
packing of 6 lines. Table 4 shows similar phenomena in higher dimensions. For example, the 
putatively best packing of 48 lines in G(9, 1) is so good that we cannot do better even when 
up through 8 lines are omitted from it. 

In the rest of this section we discuss some of the entries in Table f| that have the largest 
symmetry groups. 

Table 4 about here 



Constructions from lattices. Let Pi,. . . ,P a be a set of mutually touching spheres in a 
d-dimensional lattice packings. If there are 2N further spheres in the packing, each of which 
touches all of Pi, ... , P a , their centers form a (d — a + l)-dimensional antipodal spherical code 
with angular separation arcsec(a + 1) (||], p. 340, Theorem 1). 

In particular, the entries labeled "1" in Table |I| are obtained from the centers of spheres 
that touch one sphere in the lattices D4, Eq and Ej. The entries labeled "2" are obtained 
from spheres that touch two spheres in the lattices -D4, D$, Eq, E-j, Eg, A9 and the nonlattice 
packing P 10 fe (cf. Q, Chap. 1, Table 1.2). We remark the set of 28 lines in G(7, 1) with angle 
arccos(l/3) obtained in this way from Eg is known to be unique: see Chapter 14, Theorem 12 
of H . This configuration of lines is derived in a different way in [^] . 

It is interesting to compare Table with the table of maximal sets of equiangular lines 



given in Lemmens and Seidel j3(| and Seidel [41]. Some arrangements appear in both tables, 
for example the set of 28 lines in G(7, 1) just mentioned. On the other hand the difference 
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between the tables can be seen in dimension 8. The maximal set of equiangular lines that exists 
in G(8, 1) has size 28, with angle arcsec 3 ( |^0| , Theorem 4.6). However, Table f| shows that 
there is a set of 28 lines (not equiangular) in G(8, 1) with minimal angle arcsec 3.000511... . In 
view of the Lemmons-Seidel result (which uses the fact that in an m-dimensional equiangular 
arrangement of N lines with N > 2m the secant must be an odd integer) our set of 28 lines 
cannot be perturbed to give an equiangular set without decreasing the minimal angle. 



Prom conference matrices. It follows from Theorem 6.3 of [32| that if a symmetric confer- 
ence matrix of order q + 1 = 2 ( mod 4) exists then there is an arrangement of q + 1 equiangular 
lines in R m , m = (q+l)/2, with d 2 = (q — l)/q, meeting the simplex bound. The corresponding 
entries are labeled "3" in the table. Our program was able to find these packings for every 
prime power q of this form below 100 except for 49 and 81. 

From diplo-simplices. The entries labeled "4" are obtained by using all the vectors of shape 
±c(n 1 , (— l) n ), where c = l/y/n(n + 1). In the notation of |J, Chap. 4, these are the minimal 
vectors in the translates of the root lattice A n by the glue vectors [1] and [—1]. They form the 
vertices of a diplo-simplex ||, and have automorphism group 2 x »4 n +i. These packings also 
meet the simplex bound. 

From codes. Let C be a binary code of length m, size M and minimal distance d which is 
closed under complementation. Writing the codewords as vectors of dbl's, we obtain a packing 
of M/2 lines in G(m, 1) with d 2 = 4d(m — d)/m? . For example, a shortened Hamming code 
of length 10 with M = 32, d = 4 gives a packing of 16 lines in G(10, 1) meeting the simplex 
bound. This construction provides a rich supply of good packings. The Nordstrom-Robinson 



code [p.7|1 , f21|l , for example, yields a packing of 128 lines in G(16, 1) with d c = 15/16. 



40 lines in G(4, 1). The 80 antipodal points are (^ +1/2 ,0), (0,^ +1 / 2 ), {a^,b£ 2u+1 ), 
(a^+\b^ 2u ), {b^,a£ 2u ), {b^ +1 ,a£ 2u+1 ), where f = e 2 ^/ 8 , = 0,...,3, a = 2' 1 / i , 
b = VI - a 2 . The group of these 80 points is the group ^(^16 +^16) • 2 of order 256 (see [|10|]). 
Here V m denotes a dihedral group of order m. 

22 lines in G(6, 1). We use the vertices of a hemi-cube (an alternative way to describe the 
16-line packing) together with the six coordinate axes. 
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63 and 64 lines in G(7, 1). The putatively best packing of 63 lines in G(7, 1) (just beyond 
the range of Table ||) has angular separation 60°, and is formed from the 126 minimal vectors 
of the lattice Ej. The automorphism group of this set of 126 points is the Weyl group W(EV), 
of order 2 10 .3 4 .5.7 (cf. §, Chap. 4, §8.3). 

The best packing found of 64 lines in G(7, 1) also has an unusually large group, of order 
2|W(£ , 6)| = 2 8 .3 4 .5. This packing can be obtained as follows. The largest gaps between the 
packing of 63 lines occur in the directions of the minimal vectors of the dual lattice Ej. Let 
u be such a vector, with u = 1/2 vq, vq G E7, vq ■ vq = 6. Taking a coordinate frame in which 
the last coordinate is in the vq direction, we obtain the minimal vectors of £7 in the form 
(w,0), w £ Eq, w ■ w = 2 (72 vectors), and (x, ±^/2/3), x £ E$, x ■ x = 4/3 (54 vectors). 
We now adjoin ±V6 and rescale by multiplying the last coordinate by ^3/4, obtaining our 
final configuration of vectors (0, ±3/v2), (w,0), (x,H/y/2), which is isomorphic to what the 
computer found. The rescaling has compressed the Ej lattice in the vq direction until the angle 
between vq and the (x, l/v2) vectors is the same as the minimal angle between the (x, l/v2) 
and (w,0) vectors. This angle, the minimal angle in the packing, is arccos( \/3/ll) = 58.5178°. 

36 lines in G(8, 1). The 72 antipodal points consist of all permutations of the two vectors 
c(7 2 , -2 7 ) and c(-7 2 ,2 6 ), where c = 1/V126. These are the minimal vectors in the translates 
of the lattice by the glue vectors [2] and [—2]. 

The configuration of 16 lines in G(5, 1) is similarly obtained from the translates of A§ by 
[1], [3] and [5]. 

39 lines in G(12, 1). Take the 13 lines of a projective plane of order 3, described by vectors 
of four l's and nine 0's, and replace two of the l's by — l's in all possible ways, obtaining 78 
vectors in E 12 with angle arcsec 4. The group is 2 x £3(3), of order 2 5 3 3 13. 

If instead we replace any odd number of the l's by —l's, we obtain a putatively optimal 
packing of 52 lines in G(13, 1), with angle arcsec 4 and group 2 x £3(3). 2 of order 2 6 3 3 13. 

As in all these examples, we were amazed that the program was able to discover such 
beautiful configurations. 
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Table 3: Best packings found of N < 50 planes in G(m, 2), m < 10. The entry gives d%. 



N \ m 




4 




5 




6 




7 




8 




9 




10 


3 


1 


.5000 


1 


.7500 


2 


.0000 


2 


.0000 


2 


.0000 


2 


.0000 


2 


.0000 


4 


1 


.3333 


1 


.6000 


1 


.7778 


1 


.8889 


2 


.0000 


2 


.0000 


2 


.0000 


5 


1 


.2500 


1 


.5000 


1 


.6667 


1 


.7854 


1 


.8750 


1 


.9375 


2 


.0000 


6 


1 


.2000 


1 


.4400 


1 


.6000 


1 


.7143 


1 


.8000 


1 


.8667 


1 


.9200 


7 


1 


.1667 


1 


.4000 


1 


.5556 


1 


.6667 


1 


.7500 


1 


.8148 


1 


.8667 


8 


1 


.1429 


1 


.3714 


1 


.5238 


1 


.6327 


1 


.7143 


1 


.7778 


1 


.8286 


9 


1 


.1231 


1 


.3500 


1 


.5000 


1 


.6071 


1 


.6875 


1 


.7500 


1 


.8000 


10 


1 


.1111 


1 


.3333 


1 


.4815 


1 


.5873 


1 


.6667 


1 


.7284 


1 


.7778 


11 


1 


.0000 


1 


.3200 


1 


.4667 


1 


.5714 


1 


.6500 


1 


.7111 


1 


.7600 


12 


1 


.0000 


1 


.3064 


1 


.4545 


1 


.5584 


1 


.6364 


1 


.6970 


1 


.7455 


13 


1 


.0000 


1 


.2942 


1 


.4444 


1 


.5476 


1 


.6250 


1 


.6852 


1 


.7333 


14 


1 


.0000 


1 


.2790 


1 


.4359 


1 


.5385 


1 


.6154 


1 


.6752 


1 


.7231 


15 


1 


.0000 


1 


.2707 


1 


.4286 


1 


.5306 


1 


.6071 


1 


.6667 


1 


.7143 


16 


1 


.0000 


1 


.2000 


1 


.4210 


1 


.5238 


1 


.6000 


1 


.6593 


1 


.7067 


17 


1 


.0000 


1 


.2000 


1 


.4127 


1 


.5179 


1 


.5937 


1 


.6528 


1 


.7000 


18 


1 


.0000 


1 


.1909 


1 


.4048 


1 


.5126 


1 


.5882 


1 


.6471 


1 


.6941 


19 





.9091 


1 


.1761 


1 


.3948 


1 


.5078 


1 


.5833 


1 


.6420 


1 


.6889 


20 





.9091 


1 


.1619 


1 


.3888 


1 


.5026 


1 


.5789 


1 


.6374 


1 


.6842 


21 





.8684 


1 


.1543 


1 


.3821 


1 


.4987 


1 


.5750 


1 


.6333 


1 


.6800 


22 





.8629 


1 


.1419 


1 


.3333 


1 


.4912 


1 


.5714 


1 


.6296 


1 


.6762 


23 





.8451 


1 


.1332 


1 


.3333 


1 


.4859 


1 


.5680 


1 


.6263 


1 


.6727 


24 





.8372 


1 


.1251 


1 


.3326 


1 


.4790 


1 


.5638 


1 


.6232 


1 


.6696 


25 





.8275 


1 


.1178 


1 


.3229 


1 


.4725 


1 


.5594 


1 


.6204 


1 


.6667 


26 





.8144 


1 


.1113 


1 


.3151 


1 


.4666 


1 


.5556 


1 


.6177 


1 


.6640 


27 





.8056 


1 


.1045 


1 


.3071 


1 


.4606 


1 


.5556 


1 


.6154 


1 


.6615 


28 





.8005 


1 


.0989 


1 


.2987 


1 


.4531 


1 


.5556 


1 


.6118 


1 


.6593 


29 





.7889 


1 


.0937 


1 


.2887 


1 


.4286 


1 


.5455 


1 


.6083 


1 


.6571 


30 





.7809 


1 


.0875 


1 


.2804 


1 


.4234 


1 


.5398 


1 


.6049 


1 


.6552 


31 





.7760 


1 


.0822 


1 


.2675 


1 


.4167 


1 


.5342 


1 


.6011 


1 


.6527 


32 





.7691 


1 


.0766 


1 


.2588 


1 


.4106 


1 


.5304 


1 


.5978 


1 


.6501 


33 





.7592 


1 


.0722 


1 


.2526 


1 


.4038 


1 


.5244 


1 


.5935 


1 


.6476 


34 





.7549 


1 


.0671 


1 


.2447 


1 


.3978 


1 


.5216 


1 


.5911 


1 


.6448 


35 





.7489 


1 


.0640 


1 


.2430 


1 


.3915 


1 


.5158 


1 


.5893 


1 


.6431 


36 





.7477 


1 


.0596 


1 


.2345 


1 


.3843 


1 


.5086 


1 


.5885 


1 


.6414 


37 





.7286 


1 


.0519 


1 


.2272 


1 


.3784 


1 


.5000 


1 


.5816 


1 


.6364 


38 





.7198 


1 


.0462 


1 


.2239 


1 


.3722 


1 


.5000 


1 


.5768 


1 


.6334 


39 





.7095 


1 


.0404 


1 


.2206 


1 


.3663 


1 


.5000 


1 


.5730 


1 


.6305 


40 





.7066 


1 


.0355 


1 


.2149 


1 


.3606 


1 


.5000 


1 


.5695 


1 


.6289 


41 





.6992 


1 


.0302 


1 


.2115 


1 


.3575 


1 


.5000 


1 


.5665 


1 


.6242 


42 





.6948 


1 


.0256 


1 


.2079 


1 


.3550 


1 


.5000 


1 


.5639 


1 


.6213 


43 





.6844 


1 


.0207 


1 


.2037 


1 


.3449 


1 


.5000 


1 


.5616 


1 


.6192 


44 





.6831 


1 


.0159 


1 


.2007 


1 


.3407 


1 


.5000 


1 


.5588 


1 


.6174 


45 





.6809 


1 


.0122 


1 


.1969 


1 


.3354 


1 


.4839 


1 


.5570 


1 


.6160 


46 





.6793 


1 


.0078 


1 


.1941 


1 


.3321 


1 


.4821 


1 


.5480 


1 


.6099 


47 





.6732 


1 


.0042 


1 


.1924 


1 


.3265 


1 


.4574 


1 


.5412 


1 


.6072 


48 





.6667 


1 


.0001 


1 


.1907 


1 


.3232 


1 


.4490 


1 


.5353 


1 


.6045 


49 





.6667 





.9960 


1 


.1873 


1 


.3209 


1 


.4430 


1 


.5304 


1 


.6019 


50 





.6667 





.9910 


1 


.1841 


1 


.3150 


1 


.4370 


1 


.5257 


1 


.6003 
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Table 4: Maximal angular separation found for N < 50 lines in G(m, 1), m < 9. 



N\m 




3 




4 




5 




6 




7 




8 




9 


3 


90. 


.0000 


90 


.0000 


90 


.0000 


90 


.0000 


90 


.0000 


90 


.0000 


90 


.0000 


4 


70. 


.5288 2 


90. 


.0000 


90. 


.0000 


90. 


.0000 


90. 


.0000 


90 


.0000 


90 


.0000 


5 


63 


,4349 6 


75. 


,5225 4 


90. 


.0000 


90 


.0000 


90. 


.0000 


90 


,0000 


90 


.0000 


6 


63 


,4349 3 


70 


.5288 2 


78 


.4630 4 


90 


.0000 


90 


,0000 


90. 


,0000 


90 


.0000 


7 


54. 


,7356 6 


67. 


.0213 


73 


.3689 


80 


.4059 4 


90. 


,0000 


90 


,0000 


90 


.0000 


8 


49. 


6399 6 


65. 


.5302 


70 


.8039 


76 


.0578 


81. 


,7868 4 


90. 


,0000 


90 


.0000 


9 


47. 


.9821 6 


64 


.2619 


70 


.5288 


73 


.8437 


78. 


,4630 


82. 


.8192 4 


90 


.0000 


10 


46. 


6746 6 


64 


.2619 


70 


.5288 2 


73 


.6935 


76. 


,3454 


79. 


.4704 


83 


,6206 4 


11 


44. 


,4031 


60 


.0000 


67 


.2543 


71 


.5651 


75. 


,0179 


77. 


.8695 


80 


.6204 


12 


41. 


,8820 6 


60 


.OOOO 1 


67 


.0213 


71 


.5651 


74. 


,1734 


76. 


.6050 


79 


.4704 


13 


39 


.8131 


55 


.4646 


65 


.7319 


70 


.5288 


73 


,8979 


76. 


.1645 


77. 


.9422 


14 


38. 


.6824 


53 


.8376 


65 


.7241 


70 


.5288 


73 


,8979 3 


75. 


.0349 


77 


.2382 


15 


38. 


1349 6 


52 


.5016 


65 


.5302 


70 


.5288 


71. 


,5678 


74. 


.3318 


76 


.5006 


16 


37 


.3774 6 


51 


.8273 


63 


,4349 s 


70 


.5288 2 


70. 


,9861 


74. 


.1005 


75. 


.9638 


17 


35. 


.2353 


50 


.8870 


61 


.2551 


68 


.1088 


70. 


,5926 


73. 


,1371 


75 


.9638 


18 


34. 


.4088 


50 


.4577 


61 


.0531 


67 


.3744 


70. 


,5527 


72. 


.7464 


75 


.9638 3 


19 


33 


.2115 


49 


.7106 


60 


.0000 


67 


.3700 


70. 


,5288 


72. 


.0756 


74 


.4577 


20 


32 


.7071 


49 


.2329 


60 


.0000 6 


67 


.0996 


70. 


,5288 


71. 


.6706 


74 


.2278 


21 


32 


.2161 


48 


.5479 


57 


.2025 


67 


.0213 


70. 


,5288 


71. 


.3521 


73 


.7518 


22 


31 


.8963 


47 


.7596 


56 


.3558 


65 


.9052 5 


70. 


,5288 


71. 


.0983 


73 


.1894 


23 


30 


,5062 


46 


.5104 


55 


.5881 


63 


.6744 


70. 


,5288 


70. 


.7720 


72 


.7488 


24 


30 


.1628 


46 


.0478 


55 


.2279 


63 


.6122 


70. 


,5288 


70. 


.6027 


72 


.6547 


25 


29. 


.2486 


44 


.9471 


54 


.8891 


62 


.4240 


70. 


,5288 


70. 


.5490 


72. 


.3124 


26 


28. 


.7126 


44 


.3536 


54 


.2116 


61 


.7377 


70. 


,5288 


70. 


.5432 


72 


.1763 


27 


28. 


.2495 


43. 


.5530 


53. 


.5402 


61. 


.4053 


70. 


,5288 


70. 


.5392 


71 


.6650 


28 


27. 


.8473 


43 


.1566 


53. 


,2602 


60. 


.5276 


70. 


,5288 2 


70. 


.5322 


71. 


.5794 


29 


27. 


.5244 


42. 


.6675 


53. 


,0180 


60. 


.1344 


66 


,7780 


70. 


,5288 


71 


,5175 


30 


26. 


.9983 


42 


.2651 


52. 


.7812 


60 


.0213 


65, 


,7563 


70. 


.5288 


71 


.5175 


31 


26. 


,4987 


42. 


.0188 


52. 


.4120 


60. 


.0000 


65, 


,1991 


70. 


.5288 


70. 


.8508 


32 


25. 


,9497 


41 


.9554 


52. 


.3389 


60 


.0000 


64, 


,7219 


70. 


.5288 


70. 


.7437 


33 


25. 


.5748 


41 


.4577 


52. 


.2465 


60 


.0000 


64. 


.6231 


69. 


3203 


70 


.6940 


34 


25. 


.2567 


40. 


.9427 


51. 


.8537 


60. 


.0000 


64. 


.6231 


69. 


.1688 


70 


.6512 


35 


24. 


.8702 


40. 


.7337 


51 


.8273 


60. 


.0000 


64. 


.6231 


69. 


.0752 


70 


.6337 


36 


24. 


.5758 


40. 


.6325 


51. 


.8273 


60. 


.0000 1 


64. 


.6231 


69. 


,0752 s 


70 


.5864 


37 


24. 


.2859 


40. 


.4486 


51 


.8273 


57. 


.6885 


62. 


.3797 


67. 


.7827 


70 


.5695 


38 


24. 


.0886 


40. 


.4419 


50. 


.3677 


57. 


.1057 


62. 


.1435 


67. 


,3835 


70 


.5571 


39 


23 


.8433 


39. 


.6797 


50. 


.0611 


56. 


.8357 


61. 


.7057 


67. 


.1387 


70 


.5443 


40 


23 


.3293 


39. 


.0236 5 


49. 


.5978 


56. 


.0495 


61. 


3792 


66. 


.3815 


70 


.5288 


41 


22. 


.9915 


38. 


.5346 


49. 


.2600 


55. 


.8202 


61. 


.1630 


65. 


.9282 


70 


.5288 


42 


22. 


,7075 


38. 


.3094 


48. 


.6946 


55. 


.6160 


60. 


.8232 


65. 


.8166 


70 


.5288 


43 


22. 


,5383 


37. 


.7833 


48. 


.4030 


55. 


.3981 


60. 


.5193 


65. 


.2885 


70 


.5288 


44 


22. 


.2012 


37. 


.3474 


48. 


.0955 


55. 


.1259 


60. 


.3623 


65. 


.2422 


70 


.5288 


45 


22. 


.0481 


37. 


.1198 


47. 


.7723 


54. 


.9980 


60. 


,2282 


64. 


,7476 


70. 


.5288 


46 


21. 


.8426 


36 


.9997 


47. 


.3753 


54 


.9858 


60. 


.1101 


64. 


.4007 


70 


.5288 


47 


21. 


.6609 


36 


.5952 


47. 


.0323 


54. 


.7356 


60. 


.0433 


64. 


.1542 


70 


.5288 


48 


21. 


.4663 


36 


.3585 


46. 


.7105 


54 


.5940 


60. 


.0116 


63. 


.8846 


70 


,5288 2 


49 


21. 


.1610 


36 


.1369 


46. 


.4345 


54. 


.5031 


60. 


.0000 


63 


.4849 


68 


.0498 


50 


20. 


.8922 


36 


.0754 


46. 


.1609 


54 


.3191 


60. 


0000 1 


63. 


.1527 


67 


,7426 



Key. 1,2: from sphere packings; 3: from conference matrices; 4: diplo-simplex; 5: described 
below; 6: described in Section 3. 
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